#!/usr/bin/env python3
# -*- coding: utf-8 -*-
"""
Created on Mon Jul 17 15:00:47 2023

@author: liqingsimac
"""

'''
import numpy as np
from scipy.optimize import least_squares

x0=np.array([245,164,192,232])
y0=np.array([442,480,281,300])
d=np.array([126.2204,120.7509,90.1854,101.4021])
fx=lambda x:np.sqrt((x0-x[0])**2+(y0-x[1])**2)-d
s=least_squares(fx, np.array([0.5,0.5]))
print(s,'\n','---------','\n',s.x)
'''


import numpy as np
import matplotlib.pyplot as plt

x0=np.array([245,164,192,232])
y0=np.array([442,480,281,300])

x=149.5089 
y=359.9848

plt.plot(x0,y0,'bo')
plt.plot(x,y,'ro')
for k in range(4):
    plt.plot([x0[k],x],[y0[k],y],'b-')







